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Abstract. Isotropic Gaussian random fields on the sphere are characterized by Karhunen- 
Loeve expansions with respect to the spherical harmonic functions and the angular power 
spectrum. The smoothness of the covarianco is connected to the decay of the angular power 
spectrum and the relation to sample Holder continuity and sample differentiability of the 
random fields is discussed. Rates of convergence of their finitely truncated Karhunen-Loeve 
expansions in terms of the covariance spectrum are established, and algorithmic aspects of 
fast sample path generation via fast Fourier transforms on the sphere are indicated. The 
relevance of the results on sample path regularity for isotropic Gaussian random fields and the 
corresponding lognormal random fields on the sphere for several models from environmental 
sciences is indicated. Finally, the stochastic heat equation on the sphere driven by additive, 
isotropic Wiener noise is considered and strong c:onvcrgence rates for spectral discretizations 
based on the spherical harmonic functions arc proven. 



1. Introduction 

Path regularity of Gaussian random fields (GRFs) on subsets of Euclidean space is well 

studied, where the spectral theory of these fields is used (see, e.g., [26, 27, 28]). However, 
the general theory of second order random fields as developed in [26, 27, 28] requires a group 
structure on the space of realizations. The (practically relevant) case of GRFs indexed by the 
sphere, which we denote by S^, (and, more generally, S^") takes a special role with regard to 
invariance under (topological) group actions (see, e.g., [19] and the references there for a lucid 
discussion), so that the general results in [26] do not apply directly. Due to the relevance of 
GRFs on in applications, in particular in environmental modeling, it is of some interest to 
develop a theory of sample path regularity, stochastic partial differential equations, and their 
numerical analysis. The contribution of some basic results with direct proofs is the purpose 
of the present paper. 

Specifically, we derive the connection between the smoothness of the covariance kernel of 
an isotropic GRF on and the decay of its angular power spectrum and characterize its 
P-a.s. sample path Holder continuity and sample differentiability. Furthermore, we construct 
isotropic Q- Wiener processes using isotropic GRFs. We solve the stochastic heat equation 
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on driven by isotropic Q- Wiener noise with a series expansion with respect to the spherical 
harmonic functions. We show that the convergence rate of the fuhy discrete approximation 
scheme given by the truncation of the series expansion depends only on the decay of the 
angular power spectrum and that it is independent of the chosen space and time discretization. 

The outline of this paper is as follows: In Section 2 wc recapitulate basic definitions of 
isotropic GRFs on and of the Karhunen-Loeve expansions in spherical harmonic functions 
of these fields from [10]. A characterization of the the decay of the angular power spectrum of 
isotropic GRFs in terms of a regularity of the covariance kernel in a scale of weighted Sobolev 
spaces on §^ is presented in Section 3. Section 4 contains a version of the Kolmogorov- 
Chentsov theorem for random fields on and therefore sample Holder continuity of random 
fields is addressed. Sufficient conditions on the angular power spectrum are presented for P-a.s. 
sample Holder continuity and differentiability of isotropic GRFs. In Section 5 we approximate 
isotropic Gaussian random fields by finite truncation of their Karhunen-Loeve expansions. 
We discuss convergence rates of these approximations in p-th. moment and in the P-a.s. sense. 
The topic of Section 6 is the introduction of the practically important case of lognormal 
random fields. These are crucial in a number of applications, in particular in meteorology 
and in climate modeling. In this section, we give analogous results to Section 4, i.e., sample 
path regularity of lognormal random fields in terms of Holder continuity and differentiability 
is addressed. Finally, isotropic Q- Wiener processes are introduced in Section 7. We consider 
the stochastic heat equation on §^ driven by an isotropic Q- Wiener process and solve the 
stochastic partial differential equation (SPDE) with spectral methods. We approximate the 
solution by truncation of the derived spectral representation and show convergence rates in 
p-th moment as well as P-almost surely. These results are illustrated by numerical examples. 



2. Isotropic Gaussian random fields on the sphere 

In this section we introduce isotropic Gaussian random fields and their properties. We 
focus especially on Karhunen-Loeve of these random fields. In doing so, we follow closely 
the introduction of Gaussian random fields in Chapter 5 in [10]. Throughout, we denote by 
{Cl,A,f) a probability space and write for the unit sphere in M^, i.e., 

§2 = {x e M■^ ||x|| = 1}, 

where || • || denotes the Euclidean norm. Let (S^,(i) be the compact metric space with the 
geodesic metric given by 

d{x, y) = arccos(x, y)^3 
for all re, ?/ G S^. We denote by ;S(§^) the Borel a-algebra of S^. 

Definition 2.1. A. A® ;B(S^)-measurable mapping T : x ^ M is called a real-valued 
random field on the unit sphere. 

The random field T is called strongly isotropic if, for all G N, xi, . . . , G S^, and for 
g G 50(3) the multivariate random variables (T(a;i), . . . , T(a;^.)) and (T(gxi), . . . ,T{gxk)) 
have the same law, where 50(3) denotes the group of rotations on 

It is called n-weakly isotropic for n > 2 if E(|r(x)|"^) < +00 for all x G and if for 
1 < A; < n, xi, . . . ,Xfe G and g G 50(3), 



E(r(xi) • • • r(xfc)) = E{T{gxi) ■ ■ ■ T{gxk)). 
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Furthermore, it is called Gaussian if for all G N, xi, . . . , G the multivariate ran- 
dom variable {T{xi), . . . ,T{xk)) is multivariate Gaussian distributed, i.e., aiT{xi) is a 
normally distributed random variable for all G M, i = 1, . . . , /c. 

In what follows, we focus on real-valued random fields. Similarly to a Gaussian random 
field (GRF for short) on W^, d G N, a GRF on §^ has the following property proven, e.g., in 
Proposition 5.10(3) in [10]: 

Proposition 2.2. Let T he a GRF on S^. Then, T is strongly isotropic if and only if T is 
2-weakly isotropic. 

A key role in our analysis and simulation of isotropic GRFs on is taken by their 
Karhunen-Locve expansions. To introduce Karhunen-Loeve expansions of isotropic GRFs 
(and the corresponding Q- Wiener processes on §^ in the formulation of SPDEs on §^ in Sec- 
tion 7), we first have to define the spherical harmonic functions on which take a crucial 
role. We recall that the Legendre polynomials (P^, i G No) are for example given by Rodrigues' 
formula (see, e.g., [21]) 

for all £ G No and ^ G [—1,1]. The Legendre polynomials define the associated Legendre 
functions {Pem, ^ G No, m = 0, . . . , ^) by 

am 

PiM := (-l)-(l -;x2)W2_p^(^) 

for £ G No, m = 0, . . . and /j, G [—1, 1]. Here and throughout, we do not separate indices 
for doubly subscripted functions and coefEcients by a comma, with the understanding that 
the reader will recognize double indices as such. With this in mind, we further introduce 
the surface spherical harmonic functions y := {Yirn-,^ G No,m = —£,...,£) as mappings 
Y^jn ■ [0, tt] X [0, 2tt) — >■ C, which are given by 

for ^ G No, m = 0, . . . , ^, and (i?, (p) G [0, vr] x [0, 2tt), and by 

'■= (-l)"'^£-m, 

for ^ G N and m = —£,..., —1. By the Peter-Weyl theorem (see, e.g.. Proposition 3.29 
in [10]), y is an orthonormal basis of L^(S^;C) which we abbreviate by L^(S^). Every real- 
valued function / in L^(S^) admits the spherical harmonics series expansion 

oo e 

f = ^2 ^2 fi-m-^^rn 

e=o m=-e 

and the coefficients satisfy (cf., e.g.. Remark 3.37 in [10]) 

ftm = {—^)^fe-m, 

i.e., / can be represented in L^(S^) by the series expansion 

oo i 
/ = XI {fioYiO + 2 X ~ ^™ firn Im y^rn)) • 

e=o m=l 
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In what follows, we set for y G 

yemiy) •■= yem{^,v), 

where y = (sin cos sin i? sin cost?), i.e., we identify (with a slight abuse of notation) 
Cartesian and angular coordinates of the point y G Furthermore we denote by a the 
Lebesgue measure on the sphere which admits the representation 

da{y) = sin-dd-ddif 

for y G , y = (sin i9 cos if, sin sin if, cos i?) . 

We define the spherical Laplacian, also called Laplace-Beltrami operator, in terms of spher- 
ical coordinates similarly to Section 3.4.3 in [10] by 

It is well-known (see, e.g.. Theorem 2.13 in [13]) that the spherical harmonic functions y are 
the eigenfunctions of A§2 with eigenvalues {—^{^ -|- 1),^ G No), i.e., 

for all £ G No, m = -£,... Furthermore it is shown in Theorem 2.42 in [13] that L^(S^) 
has the direct sum decomposition 

oo 

l2(S2)=07^,(S2)^ 

where the spaces {T-L(,,l G Nq) are spanned by spherical harmonic functions 

^^(S^) = span{Y^^, m = 

i.e., H^(S^) denotes the space of eigenfunctions of that correspond to the eigenvalue 
+ 1) for I G No- 

The significance of the spherical harmonic functions lies in the fact that every 2-weakly 
isotropic random field admits a convergent Karhunen-Loeve expansion. The following result, 
which is proven in Theorem 5.13 in [10], makes this precise. 

Theorem 2.3. Let T he a 2-weakly isotropic random field on S^, then the following statements 
hold true: 

(1) T satisfies F-almost surely 

/ T{xf da (x) < +00. 

(2) T admits a Karhunen-Loeve expansion 

oo I 

(1) ^ X/ ^ '^e.mXim 

with 

aem= T{y)Yemiy)da{y), 

fori G No and m G {-^, 
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(3) The series expansion (1) converges in L^(n x §^;M), i.e., 

ah t 
{T{y) - V V aemYem{y)f da{y)) = 0. 



=0 m=-e 

■2/ 



(4) The series expansion (1) converges in L (ri;M) for all x G S , i.e., /or a/i x G § 

L ^ 

lim E((r(x) - V V afml£m(x))2) = 0. 
e=o m=-i 

This result implies that every 2-weakly isotropic random field is an element of L^(f2; L^(S^)). 
For the efficient computational simulation of 2-weakly isotropic Gaussian random fields, which 
we will call in the following just isotropic Gaussian random fields, we will exploit special 
properties of the random coefficients A = (a^^,^ G No,m = —i, . . . It turns out that the 
properties are similar to those of invariant GRFs on the torus with Fourier series expansions 
(see, e.g., [8]). First of all, we have by Remark 6.4, Proposition 6.6, and Equation (6.6) in [10] 
the following lemma: 

Lemma 2.4. Let T be a strongly isotropic random field on S^. The elements of the sequence A 
are, except for oqq, centered random variables, i.e., ¥.{aim) = for all £ & N and m = 
Furthermore they are complex-valued random variables that satisfy 



for £i, ^2 £ N and rui = —ii, ... i = 1,2, where 5nm = 1 if n = m and zero else. For the 
first element qqq, it holds that 

E(aooa^) = (^0 + E(aoo)^)<5o^<5om- 

The sequence of nonnegative real numbers {Ai,i G No) is called the angular power spectrum 
of T. 

The random variables aim CLiT'd a^-m satisfy for i EN and m = !,...,£ that 

In the case of interest in this manuscript that T is an isotropic GRF, Theorem 6.12 in [10] 
implies that A+ = {aim,£ £ No,Tn = 0, . . . , i) is a sequence of independent, complex-valued, 
Gaussian random variables. By Proposition 6.8 in [10], the elements of A+ for m ^ satisfy 
that Rea^^ and Ima^^ are symmetric random variables that are equal in law, uncorrelated, 
i.e., E(Rea£^ Ima^^) = 0, and that have variance 

E{{Re aemf) = E((Ima^^)') = Ae/2 . 

By Lemma 2.4 the elements of A that are not in A+ can be calculated by the formula 

Rea^^ = {-l)'^Reae-m, Ima^^ = (-1)™+^ Ima^_^, 

for £ G N and m = —I, . . . , —1. Furthermore we deduce from Proposition 6.11, Proposition 6.6, 
and Equation (6.12) in [10] and Lemma 2.4 in this manuscript that Rea^o is M{0,Ai) dis- 
tributed, i.e., it is normally distributed with mean zero and variance A^, and Ima^o = for 
i eN and that Reaoo is AA(E(r)2i/7r, Aq) distributed while Imaoo = 0. 
So, in conclusion, we have the following corollary. 
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Corollary 2.5. Let T be a 2-weakly isotropic Gaussian random field on S^. Then T admits 
the Karhunen-Loeve expansion 

oo e 

r = ^ ^ aemytm, 

e=o m=-e 

where (1^^,^ G No,m = —£, ...,£) is the sequence of spherical harmonic functions and the 
sequence A = {aim, ^ € No, = —i, ■ ■ ■ J) is a sequence of complex-valued, centered, Gaussian 
random variables with the following properties: 

(1) A+ = {aimj^ £ No,?77, = 0, . . . ,£) is a sequence of independent, complex-valued Gauss- 
ian random variables. 

(2) The elements o/ A+ with m > satisfy Kea^m o,nd Im a£m o-f^ independent and 
M{<d,A(,/2) distributed. 

(3) The elements o/A+ with m = are real-valued and the elements Rea^o are A/^(0, A^) 
distributed for i G N while Reaoo is J\f {^{T)2^ , Aq) distributed. 

(4) The elements of A with m < are deduced from those of A^ by the formulae 

Reaem = (-l)'"R.ea^_^, Ima^^ = (-1)'""'"^ Ima^-^. 

3. Decay of the angular power spectrum 

The error in a K-term truncation of the Karhunen-Loeve expansion of an isotropic GRF T 
on is closely related to the decay of the angular power spectrum of T. As we show next, the 

decay of the angular power spectrum is characterized by the behavior of the covariance kernel 
function that characterizes the isotropic GRF T. Often the kernel function kT is prescribed 
in applications. 

To specify this relation, we start with the definition of the kernel kx of the covariance of 
an isotropic Gaussian random field with prescribed angular power spectrum {A(,i G No). It 
is given for x, y G by the formula 

°° ^ °° 2£ + 1 

kT{x,y) := K{T{x)T{y)) = J2^iJ2 Yem{x)Yem{y) = J2Ai^—Pe{{x,y)^3). 

e=o m=-e e=o 

We observe that the covariance kernel kx just depends on the inner product resp. the (spheri- 
cal) distance. Accordingly, we denote by /c : [0, tt] — )• M the kernel as a function of the distance 
r = d{x,y), i.e., 

, , , . 2£ + 1 , , 
:= 2^Ae— Pi[cosr), 

for r G [0, tt]. A third way to look at the kernel is in terms of the inner product {x,y)^3. 
Therefore, we define fc/ : [— 1, 1] ^ M by 

ki{fi) := A;(arccos//) 

for all G [— 1, 1]. This implies overall for x,y G §^ that 

kT{x,y) = k{d{x,y)) = ki{{x,y)T^3). 

In the following Proposition, we give a relation between the weak weighted differentiability 
of the kernel kj and the decay of the corresponding angular power spectrum in terms of 
summability. 
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Proposition 3.1. For every n G No, it holds that the sequence > n) is in £^(No) 

if and only if (1 - iJ.^)"'/^^ki{iJ,), fi G (-1, 1) is in 1), i.e., 



22^ + 1 



52n 



(47r) 



< +00 



if and only if 



{l-n'^)'^dn < +00. 



Proof. Let us observe first that by definition 



/: 



^ , 2£ + 1 2f + 1 



47r 



By {P^"'^\£ G No) we denote the Jacobi polynomials given, e.g., by Rodrigues' formula 
for £ eNq, a,l3 > -1, and /x G [-1, 1]. They satisfy that 



Since Legendre polynomials are particular instances of Jacobi polynomials for a = /3 = 0, we 
conclude by recursion that 

for every n < £. This implies that 



£( 



5" 



5// 



2"f ! 

2 (^ + n)! 



" 2^ + l(^-n)!' 

where the last equation follows from the orthogonality of the Jacobi polynomials (see, e.g., [21]) 
and 



/•I / , . s 2 o'^n+l 



£\£\ 



2£+l{£-n)\{£ + n)\' 
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In conclusion we have shown that 



d 



in 



2 oo 



since for n > ^ the n-th derivative of vanishes. To finish the proof it remains to show that 
for n < ^ there exist constants c\ (n) and C2 (n) such that 

ci(nK2"<|^<C2(nK2n_ 
This follows from Stirling's inequalities 

by writing 

{l-nY'-'^ (1 - n/£)^(W^) 

and by using the properties of the exponential function. □ 

Proposition 3.1 provides a necessary and sufficient criterion for the weighted 2-summability 
of the angular power spectrum {A^^i G No) in terms of weighted square integrability of the 
n-th weak derivatives of ki with respect to the weig ht function (l-^^^n 

can be formalized 

in the framework of weighted Sobolcv spaces, which we will cover in what follows. 

For n G No, let H^{—1, 1) C 1, 1) denote the standard Sobolev spaces. We define the 
function spaces 1, 1) as the closures of 1, 1) with respect to the weighted norms 

II • l|y"(-i,i) given by 

n 

II l|2 I |2 

ll^lly"(-i,i) •= 2^ 
j=o 

where for j G No the seminorm | • |yi(_i is defined by 



Hvj{-i,i) , ^ 



Qj 2 



With this definition, 1, l),n G No) is a decreasing scale of separable Hilbert spaces, i.e., 

L2(-1, 1) = 1) D V\-l, 1) D • • • D 1) D • • • 

By Ehrling's lemma, the norm of V'^{—1, 1) is equivalent to the first and the last element of 
the sum, i.e., 

II I|2 II ||2 I |2 

iPlly^C-i,!) — II'"IIl2(-i,i) + Mv^{-i,i) 

for all u G 1, 1). We will in the sequel not distinguish between these norms by a separate 

notation. 

In what follows we are deriving further equivalent norms of V"'{—1, 1) similarly to Propo- 
sition 3.1. Therefore let us first observe that any u G L'^{—1,1) can be expanded in the 
L2(— 1, 1) convergent Fourier-Legendre series 
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with 

for all leNo. Setting 
we obtain that 



"1 

U£ := I u{x)Pi{x)dx 



At := 2-KU(,, 



i.e., M is a valid kernel kj. In conclusion we have proven the following version of Proposi- 
tion 3.1. 

Proposition 3.2. Let u G L^(— 1, 1) and n G No 6e given. Then u G 1, 1) if and only if 

1=0 

i.e., 

i=0 ^ 

is an equivalent norm in V'^{—1, 1). 

In other words we have just shown an isomorphism between the spaces V"'{—1, 1) and the 
weighted sequence spaces := + £^'')J, G No)), where (^(1 + £^''),e G No) 

denotes the sequence of weights. Our goal is to extend this isomorphism to spaces V^{—1, 1) 
with 77 ^ Nq. We define for n < 77 < n + 1 the interpolation space 1, 1) with the real 

method of interpolation in the sense of [23] by 

1) = (y"(-i, i),v^+\-i, i))^_^ ^ 

equipped with the norms || • ||\/77(-i,i) given by 

dt 



POO 

\ufvr,(^-,,i)= t~''^^-^^\K{t,u 



where the X-functional is defined by 

K{t, u) = ^^nf ^ + 

for t > 0. 

The definition of the interpolation spaces £^ for r/ ^ No is done similarly. The interpolation 
property of the spaces (see, e.g., step 4 in proof of Theorem 1.3.3 in [23] or Proposition 2.4.1 
in [22]) implies that the spaces V^{—1, 1) and are also isomorphic. So the only thing that 
remains to show to extend Proposition 3.2 to fractional differentiation orders is the equivalence 
of the norm induced by the real method of interpolation and the sequence of weights with 
parameter 77 instead of n. This is done with the interpolation theorem of Stein- Weiss in the 
following proposition. 
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Proposition 3.3. Let u G L^(— 1, 1) and rj G M+ be given. Then u G V^{—1, 1) if and only if 

^ 2/-!- 1 

e=o ^ 

i.e., 

„ 2 92^+1, 



\\u 

is an equivalent norm in y(— 1, 1). 

Proof. The proposition is already proven for 77 G No in Proposition 3.2. So let n < r/ < n+1 for 
some n G No be given and set 9 := rj — n. Applying the interpolation theorem of Stein— Weiss 
(see, e.g., Theorem 5.4.1 in [1]), we get that the weights of £^ are given by 

It remains to show that this is equivalent to ^^^(1 But this follows immediately with 

the observation that the function x^, p G (0, 1) is concave on M_|_ and satisfies (x + yY > 
2P-^{xP + yP). □ 

For more details on this we refer to Appendix A.In conclusion we have shown in this 
section that the decay of the angular power spectrum in terms of summability is equivalent 
to the smoothness of the corresponding kernel kj in terms of elements of (a scale of) weighted 
Sobolev spaces. 

4. Sample Holder continuity and differentiability 

So far, our analysis of GRFs via the Karhunen-Loeve expansion in Section 2 focused on 
mean square properties. In this section we consider sample properties of isotropic GRFs in- 
troduced in Section 2. Specifically, we are interested how the P-almost sure Holder continuity 
of isotropic GRFs depends on the decay of the angular power spectrum {A£,£ G No) which is 
one possible characterization of isotropic GRFs on by Theorem 2.3 and Lemma 2.4. 

The following lemma relates the decay of the angular power spectrum to the Holder con- 
tinuity of the kernel k at zero. 

Lemma 4.1. Let {A^, £ G No) be the angular power spectrum of an isotropic GRF on §^ with 

00 

e=o 

for some (5 G [0,2]. Then the corresponding kernel function k satisfies that there exists a 
constant Cp such that for all r G [0, vr] 

m)-k{r)\<Cpr^ 

Proof. We observe that P^(l) = 1 for all £ G No and that the derivative of Pi{x) is bounded 
by ^^(l). Therefore 



J X 



<|l-xl^ 



oo 



\m - kir)\ < ^^,?i±l|l - P,(cosr)| 
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Furthermore we have that 

|l-P^(x)| <2. 

This imphes by interpolation that 

\i-Pi{x)\ < (^\i-xf-^^^y 2^--^ <2\i-xme{£ + i)y 

for all 7 G [0, 1]. Using this estimate we obtain that 

Air 

oo 

< (27r)-^|l - cosrP ^ Ae{2£ + !){£{£ + 1))^ 

where the scries converges if '^'^^q A^i'^'^'^^ , which holds by the made assumptions for all 
7 < /3/2. Finally we observe that 

r r 2 

|1 — cosr|= / sinxcLx <rsinr = r / cosxdx<r -1, 

Jo Jo 
which implies overall with the choice /3 = 27 that 

|A:(0) - A:(r)| < C^r^, 

where 

00 

Cfi := (27r)-i ^^(2^ + + l))^^'" 

(.=0 

This finishes the proof of the lemma. □ 

Lemma 4.1 asserts Holder continuity of k{r) near r = in terms of a criterion on the 
angular power spectrum of the isotropic GRF T, while we provided i'^ criteria in Section 3. 
To relate these criteria we first observe that for e > 

00 / 00 \ ^/'^ 

Aii^+^ < C(l + e) V2 Y Aje'+^^+' 
e=o \e=o J 

by the Cauchy-Schwarz inequality, where ^ denotes the Riemann zeta function. This implies 
with Proposition 3.3 that 



i=o 

if the kernel kj is in V^{—1, 1) for some rj > P + 1. 

Our next step is to give bounds on moments of T{x) — T{y) for x, y G in terms of the 
geodesic distance d{x,y). We prove the lemma by expressing the moments in terms of the 
kernel k and by an application of the preceding lemma. 

Lemma 4.2. Let T be an isotropic Gaussian random field on §^ with angular power spectrum 
{A£,i G No). // the angular power spectrum satisfies that 



J2Aei'+^<+oo 
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for some fi G [0, 2], then for all p there exists a constant Cp^p such that 

E{\T{x)-T{y)fP)<Cp,pd{x,yfP 

for all X, y G S^. 

Proof. First note that T{x) — T{y) is a centered Gaussian random variable. Furthermore, if 
X is a M{0, (T^)-distributed random variable, then 

K{\XfP) = K{\aYfP) = {a^)PE{\YfP) = K{X'^)Pc2p 

for p G N, where y is a standard normally distributed random variable and C2p denotes the 
2p-th moment of Y. We also observe that E(|T(a:;) — r(y)p) can be expressed in terms of k 
since 

E(|r(x) - r(y)|2) = E(r(,x)2) - 2E(r(z)T(y)) + E{T{yf) 
= krix, x) - 2kT{x, y) + kriy, y) 
= 2{k{0)-k{d{x,y))). 
Combining the two previous observations, we conclude that 

E(|r(x) - Tiy)\'P) = C2,E{\T{x) - T{y)\y 
= 2c2pik{0)-k{d{x,y)))P 
<2c2pCld{x,yfP, 

where we applied Lemma 4.1 in the last step. This finishes the proof of the lemma. □ 

The following result is a version of the Kolmogorov-Chentsov theorem for random fields 
with domain S^. Note that in this result the fields do not have to be Gaussian or isotropic. 

Theorem 4.3 (Kolmogorov-Chentsov theorem). LetT he a random field on the sphere that 
satisfies for some p > and some e G (0, 1] that there exists a constant C such that 

n\nx)-T{y)\P)<Cd{x,yf^'P 

for all x,y G S^. Then there exists a continuous modification of T that is locally Holder 

continuous with exponent 7 for all 7 G (0, e) . 

Proof. Let us first construct six charts {Ui,i = 1, ... ,6) that cover the sphere by taking the 
six possible hemispheres given by the coordinate system such that the boundary is a circle 
of radius r with r G (-^2/3, 1), i.e., wc take a bit less than the complete hemispheres but 
enough to cover the whole sphere. Let the coordinate maps {ipi,i = 1, . . . , 6) be given by the 
projection onto the plane that divides the hemispheres, i.e., if U is contained in the northern 
hemisphere than the corresponding coordinate map (p is given by ip{{xi,X2,X3)) := {xi,X2) 
for X = {xi,X2,X3) G U and maps onto the disc {x G M^, ||x||ir2 < r}. 

For a given chart (U,(p), we have to show that the Euclidean norm in is equivalent to 
the metric on S^, i.e., that there exist constants Ci, C2 > such that for all x,y e U 

Ci\\<f{x) - (/7(y)||]R2 < d{x,y) < C2\\f{x) - <y5(?/)||]R2 

or equivalently that 

^ ^ arccos((a;,y)R3) ^ ^ 

^ ~ Wfi^) - v'(y)llE2 ~ ^' 

We show this estimate for U contained in the northern hemisphere. The calculations for the 
other five charts are similar and the bounds are the same due to symmetry. 
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One first calculates that 



and shows that 



{x,y)R3 = 1 - ^(||</:'(a;) - ip{y)\\l^2 + |x3 - yaf) 



2r^ 



This implies that we can bound the quotient of interest from above and below by 
arccos(l - - ^ arccos((.T, y)R3) 



\\ip{x) - ^{y)\\^2 \\^p{x) - ip{y)\\^2 

arccos(l -{\ + y^)\W{x) - if{y)\[{ 



< 



\\ip{x) - ip{y)\\^2 

since arccos is a monotonically decreasing function. Let us define / : [0, 2r) — )■ M by 

arccos(l — aa^) 



/(«) ■■- 



a 

for a G (0, 2r), where a = 1/2, 1/2 + r^/(l — r^). Then one shows with standard methods 
from real analysis that / is well-defined on [0, 2r) and monotonically increasing which leads 
with the observation that /(O) = \/2a by I'Hopital's rule to the conclusion that 

_ ^ ^ arccos((x,i/)R3) ^ arccosl ^^'+j^^-^ ) _ ^ 

- \\(p{x) - (^(y)||M2 - 2r 

and finishes the proof of the equivalence of geodetic and Euclidean distances on the sphere 
and in the charts. 

For a,b E (p{U) it holds for the random field on the chart by the made assumptions and 
the equivalence of the distances that 

Ei\T{^-\a)) - T{<f-\b))n < Cd{<p-\a),^-\b)f/P+' < C ■ Cl'^^'\\a - 6f /^'+^ 

Since ^p{U) is a domain in R^, wc obtain by the Kolmogorov-Chentsov theorem for domains 
(see Theorem 2.1 in [11] or Theorem 4.5 in [16]) that there exists a continuous modification T\o 
ip~'^ that is locally Holder continuous with exponent 7 for all 7 G (0, e) and so is Ti on U due 
to the smoothness of the coordinate map. 

With the same proof we obtain continuous modifications (Tj,i = 1,...,6) on all charts 
{Ui, i = 1, . . . , 6). We glue these together with a smooth partition of unity {pi,i = 1, . . . , 6) 
on which is subordinate to the open covering, (see, e.g.. Theorem 1.73 in [9]), by 

6 
i=l 

for all x G where Ti{x) = for x ^ J7j. Then T is a continuous modification of T that is 
locally Holder continuous with the same exponent 7 for all 7 G (0, e) due to the smoothness 
of the partitions of unity. This finishes the proof of the theorem. □ 

With the made observations, we are now prepared to prove one of the main results of this 

section which states that if the angular power spectrum of an isotropic Gaussian random field 
is summable with weights than there exists a continuous modification which is Holder 

continuous of exponent 7 for all 7 < /?/2. 
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Theorem 4.4. Let T he an isotropic Gaussian random field on with angular power spec- 
trum {A£,£ G No). If the angular power spectrum satisfies that 

oo 

J^A^'^^ <+oo 
e=o 

for some /3 G [0, 2], then there exists a continuous modification ofT that is Holder continuous 
of exponent 7 for all ^ < 13/2. 

Proof. The claim follows by the application of the previous results in the following way: It 
holds by Lemma 4.2 that for all p G N and x, y G the random field satisfies 

Ei\T{x) - Tiy)fn < C0,pd{x,y)^P = C0,pd{x,yf^W2-Vp)2p_ 

Theorem 4.3 implies that there exists a continuous modification that is locally Holder contin- 
uous of exponent 7 for all 7 < /3/2 — 1/p for any p G N, i.e., of exponent 7 for all 7 < /3/2. □ 

Just as an example let us calculate the parameters of P-almost sure Holder continuity for 
the two choices of a that we simulate in the following sections. For a = 3 we get /3 < 1 which 
implies 7 < 1/2 in Theorem 4.4 and a = 5 implies /3 = 2 and therefore 7 < 1. 

Furthermore as second main result of this section we are interested in the assumptions on 
the angular power spectrum that imply the existence of differentiable modifications of isotropic 
GRFs. In particular in the context of approximate, numerical solutions of partial differential 
equations, regularity properties of samples are essential for the derivation of convergence rates 
for, e.g.. Finite Element or Finite Difference discretizations. 

Theorem 4.5. Let T he an isotropic Gaussian random field on with angular power spec- 
trum {A£,i G No) and expectation zero. If the angular power spectrum satisfies that 

00 
e=o 

for some (3 > 0, then there exists a modification ofT that is k-times continuously differentiahle 

for all nonnegative integers k < j3/2 — 1. 

Proof. Let us first observe that the made assumptions imply that T has a continuous modifica- 
tion by Theorem 4.4. Without loss of generality let T already be the continuous modification, 
which is an isotropic Gaussian random field with the same parameters and has a Karhunen- 
Loeve expansion with the same parameters by Corollary 2.5. 

We observe that the norms in the Sobolev spaces H'^{§^) C L2(§2) for > are equivalent 
to the norm || • ||a^ given by 

00 £ 

iHii:=hiii.(s^) + ii(-A)''/^«iii.(s.) = ^ ^ \uemni+m+m 

e=o m=-e 

for u G H'^{S'^) with expansion 

00 e 
e=o m=-e 
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Furthermore it holds that H^{S'^) is continuously embedded into C*^(S^) for all A; < 77 — 1 by 
Theorem 6.5 in [25]. Now T satisfies that 

/ OD e \ 



00 

E 

e=o \ m=l 



= ^^,(^ + l)(l + (^(£ + l)n 
< +00 

for all ry < /3/2 by our assumptions, where we applied Corollary 2.5. Together with a Cauchy 
sequence argument this implies that the series expansion of T converges in L?'{Vt] H^{S'^)) and 
therefore in probability. By the Ito-Nisio theorem [4] the series converges P-almost surely to 
the same limit. Therefore T G iJ^(S^) P-almost surely and by the above Sobolev embedding 
and the continuity of the random field, T G C'^(S^) for all A; < 13/2 — 1 P-almost surely. 
Setting T to zero on the remaining P-null set of Jl, we obtain a modification of T in C'^(S^) 
for all k < 13/2 — 1. This finishes the proof of the theorem. □ 

5. Approximation of isotropic Gaussian random fields 

Let us approximate and simulate isotropic Gaussian random fields in this section, where 
we use the properties of the random fields that were introduced in Section 2. In what follows, 
we consider centered random fields without loss of generality. It is clear by Corollary 2.5 that 
we can transform the centered, isotropic random field into a field with nonzero expectation 
by adding the expectation, which is a constant according to Lemma 2.4. To prepare the 
presentation of the approximation of isotropic GRFs on S^, we rewrite its series expansions, 
where we use the properties of the spherical harmonic functions and the structure of real- 
valued random fields. 

Lemma 5.1. Let T be a strongly isotropic, centered, Gaussian random field. For £ E N, 
m = 1, . . . and 1? G [0, tt], set 



Then, for y = (sin •& cos (p, sin -d sin ip, cos 1?), 

00 e 
T{y) = Y ^eXjoLeoit)) + V^e E ^em{'&){Xj^ cos(m(^) + sm{m^)) 
e=o m=l 

in law, where £ G = 0, ...,£) is a sequence of independent, real-valued, 

standard normxdly distributed random, variables and X|q = for ^ G Nq. 

Proof. By Corollary 2.5, T can be represented in the (mean-square convergent) Karhunen- 
Loeve expansion 

00 t 

r = E E ^irnXlm- 
i=Q m=-e 
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This sum can be rewritten to 



oo I 

e=o m=l 

oo e 

e=o m=l 

oo e 

^=0 m=l 

oo ^ 

= + 2Re(a£^y£^)) 




^"=0 m=l 

by Lemma 2.4 and the properties of the spherical harmonic functions. We observe that 

Yem{.'&,(p) = Lim{-&)e''^'^ = Lirn{'d){cos{m^) + ism{rrnp)) 

for if) G [0, tt] X [0, 27r) and therefore, by the properties of complex numbers, that 

Re(a^^y^j„(i?, = L^OT(??)(Rea^OTCOs(m9?) - Ima^^sin(m(^)). 

Let {{Xj^,Xf^),£ e No,m = 0, . . . , i*) be a sequence of independent, real-valued, standard 
normally distributed random variables, then 

2 ""tm, ti/t y 2 

in law for £ G N and m = 1, . . . ,^ by Corollary 2.5. Furthermore the corollary implies that 

for £ G Nq. The insertion these observations into the Karhunen-Loeve expansion of T com- 
pletes the proof. □ 

For a given sequence {{Xj^, Xf^),£ G No, m = 0, . . . ,i) as specified in Lemma 5.1, set 

oo t 

T{y) := J2 ^/AeXjoLeo{l^) + V^e J2 ^irnmxL cosimip) + Xj^ sin(m(^)). 

(.=0 m=l 

In what follows, we truncate the series expansion in order to implement it and prove its 
convergence. For k G N, we set 

K t 

T-{y) := ^^A^X}^L,^{'^) + ^ L,„(7?)(X|„ cos(m<^) + X^^ sin(m(^)), 

i=Q m=l 

where y = (sin i? cos sin i? sin cos??) and {'d,ip) G [0,7r] x [0, 27r). 

Proposition 5.2. Let the angular power spectrum {A£,i G Nq) decay algebraically with order 
a > 2, i.e., there exist constants C > and Iq E N such that A( < C ■ 1^°' for all i > Iq. 
Then the series of approximate random fields {T'^,k G M) converges to the random field T 
in L^(f2; L^(S^)) and the truncation error is bounded by 

\\T - r'*||L2(n;L2(§2)) < C ■ K~(°~2)/2 
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for K> £o, where 

\a-2 a-lJ 

Proof. Since {{Xjj^,X^^),£ G No,m = 0, ... is a sequence of independent, standard nor- 
mally distributed random variables, the error is equal to 

1 1 1 1 

llL2(n;L2(s2)) 



J2 UE((Xjo)^)||F^o|li2(s.) 



e=K+i 



+ ^(E((X|^)2)||Rey,„||i.(g.) +E((X|j2)||imy^^||2^ 



(§2) 



m=l 



= Yl (^^ll^^olli2(§2) + 2^ ^(||ReY£^|||2(§2) + ||Imyg^|||2(§2))). 

e=K+l m=l 

We observe that ||y^o|li2(s2) = 1 and || Rey^^|||2(§2) + || Iml£^|||2(§2) = 1 for £ e No and 
m = 1, . . . ,i. Therefore, the sum simplifies to 

oo 

11^- ^"^IIl2(Q;L2(§2)) = {2£ + 1)A(, 

e=K+i 

which is bounded by 

oo oo 

e=K+i e=K+i 

due to the assumed properties of the angular power spectrum. We rewrite the sum and bound 
it by the corresponding integral which leads to 

oo oo 

J2 (2r(°-i) + r") = ^(2(£ + + {£ + 

e=K+i e=i 

POO 

< / (2(x + + (x + k)-") c/x 

Jo 

\a — 2 a — I / 

This finishes the proof since is bounded by 1. □ 

In an implementation in MATLAB we verified the theoretical results. We took as "exact" 
solution the random fields with /c = 2^ since for larger k the elements of the angular power 
spectrum A£ and therefore the increments were so small that MATLAB failed to calculate 
the series expansion. Instead of the L^(§^) error in space, we used the maximum over all grid 
points which is a stronger error. In Figure 1 the results and the theoretical convergence rates 
are shown for a = 3, 5. One observes that the simulation results match the theoretical results 
in Proposition 5.2. 

Since we discussed P-almost sure Holder continuity in Section 4, we are also interested in 
P-almost sure convergence rates of the approximate random fields {T'^,k G N). Therefore, 
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(a) Angular power spectrum with parameter a = 3. (b) Angular power spectrum with parameter a = 5. 



Figure 1. error of the approximation of Gaussian random fields with 
different angular power spectrum and 1000 Monte Carlo samples. 



we include the following result on L^(S^)) convergence since we need it for optimal 

pathwise convergence rates of the approximate random fields (T*^, k G N). 

Lemma 5.3. Let the angular power spectrum {A£,i G No) decay algebraically with order 
a > 2, i.e., there exist constants C > and £o £ N such that < C ■ for all i > Iq. 
Then the series of approximate random fields {T'^,k G N) converges to the random field T 
in U\Q.; L^(S^)) for any finite p > 0, and the truncation error is hounded by 

\\T - 7"^||LP(n;L2(§2)) < Cp ■ 

for K > io, where Cp is a constant that depends on p, C , and a. 

Proof. For p < 2 the result follows with Proposition 5.2 and Holder's inequality. Therefore 
let us consider p > 2 now. We prove the claim for p = 2m, m G N. For all other p £ M+, 
the result follows again by Holder's inequality. So let m G N, then Corollary 2.17 in [2] states 
that there exists a constant Cm such that 

\\rp rpK\\2m ^ IIT^ rpK\\2m 

Applying Proposition 5.2 we conclude that 
where 



\a — z a — 

which finishes the proof. □ 

We have just shown that the convergence rate does not depend on p. This is necessary 
to get up to an epsilon the same sample convergence rates as in the p-th moment by the 
Borel-Cantelli lemma, which we show in what follows. 
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(a) Angular power spectrum with parameter a = 3. (b) Angular power spectrum with parameter a = 5. 



Figure 2. Error of the approximation of a sample of Gaussian random fields 
with different angular power spectrum. 

Corollary 5.4. Let the angular power spectrum {Ai,l G Nq) decay algebraically with order 
a > 2, i.e., there exist constants C > and Iq £ N such that < C ■ for all i > Iq. 
Then the series of approximate random fields {T'^,k £ N) converges to the random field T 
¥-almost surely and the truncation error is bounded by 

||r-r'^||i2(§2) < K-'^, F-a.s. 

for all(3< {a -2)12. 

Proof. Let /3 < (a — 2)/2, then the Chebyshev inequality and Lemma 5.3 imply that 
P(||r - T^Wmn^) > < K^PE{\\T - T%,^^,^) < C'Pk(^-("-2)/2)p_ 

For all p > {{a - 2)/2 - the series 

oo 

^k(^-("-2)/2)p<+oo 
ft=l 

converges and therefore the Borel-Cantelli lemma implies the claim. □ 

In Figure 2, we show the corresponding error plots to Figure 1 but instead of a Monte 
Carlo simulation of the approximate L^(0; L^(S^)) error we plotted the error of one sample. 
The convergence results coincide with the theoretical results in Corollary 5.4. 

To give the reader an idea of the structure of the Gaussian random fields in dependence of 
the decay of the angular power spectrum, we include two samples in Figure 3. Here we chose 
Ai = {i+ 1)-° for £ E No and a = 3, 5. Therefore, Ai < for all £ > I, which meets the 
assumptions of Proposition 5.2. We truncated the series at k = 100, since larger k do not 
affect the pictures, but the numerical accuracy suffers due to multiplication and addition of 
very small numbers. 

We remark that similarly to fast Fourier transforms, there exist fast transforms for spherical 
harmonic functions, see, e.g., [12] and the set of C routines SpharmonicKit explained in [3]. 
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(a) Angular power spectrum with parameter a = 3. (b) Angular power spectrum with parameter a = 5. 

Figure 3. Two samples of isotropic Gaussian random fields with different 
angular power spectrum and truncation at k = 100. 

These allow to simulate isotropic Gaussian random fields with the suggested approximations 
efficiently also for large choices of k. 

6. LOGNORMAL ISOTROPIC GAUSSIAN RANDOM FIELDS 

In this section we consider lognormal random fields on S^, i.e., if T is an isotropic Gaussian 
random field on then we are interested in exp(r) given by exp(T(x)) for all x G S^. These 
random fields are especially of interest when modeling Saharan dust particles (see, e.g., [15]), 
feldspar particles (cf., [24]), and ice crystals (cf., [14]). We show in the following that the 
sample regularity of a lognormal random field is the same as that of the underlying Gaussian 
random field. This is done by first proving regularity in LP{fl;M) and then applying the 
Kolmogorov-Chentsov theorem similarly to Section 4. 

Lemma 6.1. Let T be an isotropic Gaussian random field on with angular power spectrum 
{Ai,£ £ No). // the angular power spectrum satisfies that < Ci~°' for all i € N, some 
a > 2, and some constant C, then for all p £ N and f3<a — 2, fi<2 there exists a 
constant Cp^p such that 

II exp(r(x)) - exp(r(y))||ip(Q.R) 

<2exp(pA;(0))C^,pd(x,2/)^/2 

for all x,y £ E'^. 

Proof. Let us first observe that for a, 6 G M it holds that 

< |a-6| max{e^e^} < |a - 6|(e'' + e^). 
This implies for x, y G that 

II exp(r(x)) - exp(r(y))||^,(^^j,) < E((exp(r(x)) + exp(r(y)))^|r(x) - T{y)\p) 

< E((exp(r(x)) +exp(r(y)))'^)'^' .E(|r(x) -r(2/)|2P)i/2^ 
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(a) Angular power spectrum with parameter a = 3. (b) Angular power spectrum with parameter a = 5. 
Figure 4. Corresponding lognormal samples to Figure 3 with k = 100. 

where we apphed Holder's inequality in the last step. By Lemma 4.2 the second term is 
bounded by 

E{\T{x) - r(y)pP)i/2 < cl^d{x,yr''/^ 
for any /3<a — 2,/3<2. The first term satisfies that 

E((exp(T(x)) +exp(r(y)))'^)'^' < 2(2p-i)/2 (^E(exp(2pr(a;))) + E(exp(2pT(y)))) 

Since T(x) and T{y) are real-valued Gaussian random variables with expectation zero and 
variance A:(0), the moment generating function is given by 

E(exp(2pr(x))) = exp(2/A;(0)), 



which implies that 

1 /2 

E((exp(r(x)) + exp(r(y)))^^) < I^'^P'^^'^I^^^ expip'^k{0)) = 2^ exp(p2fc(0)). 

Therefore we overall conclude that 

||exp(r(x)) -exp(r(y))|Up(j,.K) < 2 exp(p fc(0))C^,pd(x, 
which finishes the proof. □ 

The lemma enables us to conclude that the lognormal random field of an isotropic Gaussian 
random field T has the same sample Holder continuity properties as T. 

Corollary 6.2. Let T be an isotropic Gaussian random field on with angular power spec- 
trum {Ai,i G No). // the angular power spectrum satisfies that Ai < C ■ l^" for all £ G N, 
some a > 2, and some constant C , then there exists a modification of exp(T) that is Holder 
continuous of exponent s for all s < [a — 2) /2, s < 1. 

Proof. The proof is the same as the one of Theorem 4.4, where we apply Lemma 6.1 instead 
of Lemma 4.2. □ 
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In Figure 4 we took the Gaussian random field samples that are shown in Figure 3 and 
plotted the deformed sphere with the corresponding lognormal radius which is done when 
modeling dust or feldspar particles resp. ice crystals. 

In Theorem 4.5 we have shown the existence of fc-times continuously differentiable modifi- 
cations of isotropic GRFs depending on the convergence of the corresponding angular power 
spectrum. The compactness of the unit sphere, the smoothness of the exponential func- 
tion, and the chain rule imply as a direct consequence that the same properties hold for the 
corresponding lognormal random fields. 

Corollary 6.3. Let T be an isotropic Gaussian random field on with angular power spec- 
trum {A£,£ G No) and expectation zero. If the angular power spectrum satisfies that 



f^Aii^+f^ <+oo 

for some /3 > 0, then there exists a modification of the corresponding lognormal random 
field exp(T) that is k-times continuously differentiable for all nonnegative integers k < P/2—1. 

7. Stochastic partial differential equations on the sphere 

In this section we consider the heat equation on the sphere with additive Q- Wiener noise as 
an example of a stochastic partial differential equation (SPDE) on S^. To discuss stochastic 
partial differential equations we first introduce Q- Wiener processes on the sphere. 

To this end let us consider Q- Wiener processes that take values in L^(§^) and that are 
isotropic in space. Then, by Lemma 5.1 and by the construction of Q- Wiener processes out 
of GRFs as was done in an abstract setting, e.g., in [2, 17], a Q- Wiener process taking values 
in L^(S^) can be characterized by the Karhunen-Loeve expansion 

oo e 
e=o m=-e 

oo e 

= ^ ./Ae(3lo{t)Yeo{y) + E (/^^-(i) ^eYUy) + PLit) Imy,^(y)) 

e=o m=l 

oo e 

= VAef5Ut)Leo{^) + yplA, ^irnm^Lit) cos(m(^) + sin(m(^)), 

where y = (siniJcosi^, sin sin cost?) and ((/^j^^, /3|„), ^ G No,m = 0, ...,£) is a sequence 
of independent, real- valued Brownian motions and /?|q = for £ G Nq and t G M+. The 
covariance operator Q is characterized similarly to the introduction in [7] by 

oo £ oo £ 

QYlM = Y Y mW{l),YLM)HiWil),Ye^)H)Yim = Y Y E(aLM(l)^Ul))^^m 

e=0 m=-l £=0 m=-l 

= AlYlm 



for L G No and M = —L, . ■ ■ ,L, i.e., the eigenvalues of Q are characterized by the angular 
power spectrum {A£,£ G No) and the eigenfunctions are the spherical harmonic functions. 
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Let US calculate ||W^(i)llL2(n;L2(S2)) for t € M+ next. It holds similarly to the proof of 
Proposition 5.2 that 

oo 

£ 

+ 2^,5^(E((/3]„(i))')||Rey,^||i2(§2)+E((/3L(i))')l|Iml£.„|li2(s.)) 

m=l 

oo 

= t^{2£ + l)Ae = tTtQ. 
This expression is finite for any finite t G M+, if 

oo 

^A^i < +00. 

e=o 

With the made introductions of Q- Wiener processes and the Laplace operator on the sphere 
in Section 2, we are now able to write down the stochastic heat equation 

(2) dX{t) = As2X{t) dt + dW{t) 

with initial condition X{0) = Xq e L"^ {Q; L"^ , where t € T = [0,r], T < +00. 
Looking for solutions in L^(S^), we rewrite Equation (2) to 

X{t) = Xo+ [ A§2X(s) ds+ [ dW{s) =Xo+ [ A§2X(s) ds + W{t) 
Jo Jo Jo 

and further, since the spherical harmonic functions y form an orthonormal basis of -L^(S^) 
and are eigenfunctions of Ag2 , we have that 

00 e 

{X{t),Yim)L2{S'2)Yim 

= XI XI ( (^0, ^m)L2 (§2)^^,71 + / iX{s),Yijn)L'2{S.'2)'^S2Yijnds + aim{t)Yim\ 

e=o m=-e 

00 e t 

X X ( (^0, i^m)L2(S2) — -^(^ + 1) / {^{s),Y£jn)j^2(g2)ds + aim m- 

e=o m=-e 



e=o m=-e 

00 



This is equivalent to solve for all £ G No and m = . . . ,i the stochastic (ordinary) differ- 
ential equation 

Y^^)i2(§2) = {Xo, y^^)^2(s2) - £{£ + 1) / {X{s), y^„)^2(g2) ds + aem{t)- 

Jo 

The variations of constants formula yields 

(X(i),y,^)^2(§2) = e-^(^+^)*(Xo,F,„)^2(s2) + re-^(^+i)(*-^)da,^(s). 

Jo 
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So overall, the solution of the stochastic heat equation (2) reads 

oo I 



1=0 m=-e •'^ 

oo i „t 

e=o m=-e 

+ V2^( re-^(^+^)(*-)d^L(.)Rey,^+ T e-^(^+^)(*-) d^L(.) Imy,^))) 

m=l -^0 

oo 

=:E^^(*)' 

e=o 

and wc choose the sequence of stochastic processes {X£,i G Nq) accordingly. Each process in 
this sequence satisfies the recursion formula 

at+h 

^ ft+h 
m=l 

Similarly to [5] we observe that by the Ito formula (see, e.g., [6]) 

(3) f e-^('+'^('-^UpUs) 

Jo 

is normally distributed with mean zero and variance (2£{i + 1))~^(1 — e~^^(^+^)*) for £ eN, 
m = 1, . . . ,i, and i = 1,2. This implies that 

rt+h 

where 

aj„:= (2£(£ + l))-i(l-e-2W)'*). 

For ^ = we have no convolution integral and therefore the distribution of the expression is 
that of (the increment of) a standard Brownian motion, i.e., cTq^ = h. 

For the simulation of paths of the solution, we have to compute the solution on a discrete 
time grid = to < ti <■■■< tn = T, n e N, on which the path of the Brownian motion resp. 
the stochastic integral (3) is known. The stochastic integral (3) at time tk, fc = 0, . . . ,n, is 
equal in law to a sum of weighted, standard normally distributed random variables 



Jo Jtj 

k—1 I 
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j=0 

where hj = tj+i - tj, j = 0, . . . .n - 1 and {X'l^^{j)J. G No,m = 0, i = 1,2, j = 

0, . . . , n — 1) is a sequence of independent, standard normally distributed random variables. 
This enables us to write down the solution of Equation (2) recursively 

e 

= e-^^^+^^^>'Xe{tk) + ^Aeaeh, (xjo{k)Yeo + i^Lik) Re Y^^ + ImYem)) 

m=l 

for all £ G No and k = 0, . . . ,n—l. Using the notation of Lemma 5.1, we rewrite the recursion 
to 

Xe{tk+i) = e-'^'+'^^'^Xi{tk) + Mk) 

m=-e j=0 
where the increments are given by 

e 

y) = \/Ae(Jeh. (xlo(j)L,o(^) + LemW {xUj) cos(m</.) + Xf^{j) sm{mip))) 

m=l 

for y = (sin 1? cos sin?? sin cos??) G and j = 0, . . . ,n — 1. We observe for later use that 

j=0 



3=0 



m=l j=0 

+ e-^(^+i) '^ae,^XUj)) sin(m^))) 



j=0 



and that 



j=0 

is a normally distributed random variable with mean zero and variance 
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This implies that we have equality in law of 

(xloL^o W + V2 ^ LimW {Xj^ cosimcp) + Xf^ sin(m(^))) , 

m=l 

where {{X}m-> -^fm) ■> ^ No,m = 0,...,^) is a sequence of independent, standard normally 
distributed random variables. 

To implement the solution, we calculate X^ exactly for finitely many ^ G No on a finite time 
and space grid. One way to discrctizc the sphere is to take an equidistant grid in ■!? G [0, tt] 
and (p G [0, 27r). Then we add the calculated X^ and get an approximate solution, i.e., we 
simulate the approximate solution X*^, k G No by 

K 

on finitely many time and space points. In what follows let us estimate the mean square error 
when truncation of the series expansion at k G N is done. 

Lemma 7.1. Let t £ T and = to < ■ ■ ■ < tn = t be a discrete time partition for n G N, 
which yields a recursive representation of the solution X of Equation (2). Furthermore, 
assume that there exist Iq E N, a > 0, and a constant C > such that the angular power 
spectrum (Ai, i G Nq) satisfies Ai < C ■ for all i > Eq. Then the error of the approximate 
solution X'^ is bounded uniformly in time and independently of the time discretization by 

\\X{t)-X-{t)\\L2^n;LHS^))<C-K-"/^ 

for all K> £o, where 

= ll^o|li2(n.i2(§2)) + + ^qry) • 

Proof. Let t G T and = to < • • • < tn = t be a partition of [0, t] for some n G N. Since 
E{ipi{j)) = for all ^ G No and j = 0, . . . , n — 1, we first observe that 

(4) 

\\X{tn) - -'^'*(*n)||L2(n.i2(§2)) 

oo £ oo n— 1 

e=K+i m=-e e=K+i j=o 



2 

L2(Q;L2(S2)) 



£ ^ e-^(^+i)*"(Xo,y,^),.2(s2)y,^ 

?=/t+l m=-e 

oo n— 1 

+ 11 E T.-~'^'^"^^-^^'''Mj) 

e=K+i j=o 



2 

L2(f2;L2(§2)) 

2 

L2(n;L2(S2))' 



We define an isotropic Gaussian random field as in Lemma 5.1 by 

oo i 

T:=Y^ ^/A^aa^ (x^L^o(^?) + ^2 J] L^^) {Xj^ cos(m<^) + X^^ sm{m<f)) 



e=o m=l 
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with angular power spectrum {A£aj^^,£ G No) and denote similarly to Section 5 by T'^ the 
truncated series expansion. Then 



oo n—1 



'.=K+1 j=0 



2 

I j rp rJ-^K 1 1 2 



L2(f2;L2(§2)) II IIL2(Q;L2(§2)). 

The angular power spectrum satisfies with the made assumptions that 

A,al = A,— 3— (1 - e-2^(^+i)*") < cr'^r^ ■ 1 = cr("+2) . 

t u„ *2^(^ + 1)^ ^ - 

With these parameters we apply Proposition 5.2 to the difference of T and T'^ which yields 

/ 2 1 \ 

The first term in the last line of (4) is bounded by 



oo e 2 



e=K+i m=-e 



L2(f2;L2(S2)) 



£=K+1 m=—£ 

Therefore, it converges faster than any polynomial, especially it can be bounded by So 
overall we bound 

- X'^{tn)\\'i2^Q.L2(§2)) < (^C ■ (^- + a + l) 

which finishes the proof of the lemma. □ 

We remark that it is not necessary that the angular power spectrum {A£,£ G No) of the 

Wiener process decays with for a > 2 but that it is sufficient to assume that a > 0. 

In an implementation in MATLAB we verified the theoretical results of Lemma 7.1. We 
took as "exact" solution the approximate solution at time T = 1 with k = 2^ since for larger 
K the elements of the angular power spectrum and therefore the increments were so small 
that MATLAB failed to calculate the series expansion. We calculated the solution in one 
time step since we have shown in Lemma 7.1 that the convergence rate is independent of the 
number of calculated time steps. Instead of the L^(§^) error in space, we used the maximum 
over all grid points which is a stronger error. In Figure 5 the results and the theoretical 
convergence rates are shown for a = 1,3,5. One observes that the simulation results match 
the theoretical results from Lemma 7.1. 

Similarly to the proof of almost sure convergence of approximations of isotropic Gaussian 
random fields in Section 5, we need a convergence result for the approximation of the 
solution of the stochastic heat equation to show pathwise convergence. This is proven in the 
following by a combination of Lemma 5.3 and Lemma 7.1. 

Lemma 7.2. Let t eT and = to < ■ ■ ■ < tn = t be a discrete time partition for n 6 N, which 
yields a recursive representation of the solution X of Equation (2). Furthermore assume that 
there exist Iq E N, a > 0, and a constant C > such that the angular power spectrum 
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(c) Angular power spectrum with parameter a = 5. 



Figure 5. !?■ error of approximation of the stochastic heat equation different 
angular power spectrum of the Q- Wiener process and 100 Monte Carlo sam- 
ples. 



(A£,i G No) satisfies A£ < C ■ i " for all i > io. Then the error of the approximate solution 
X'^ is hounded uniformly in time and independently of the time discretization by 

\\X{t) - -^''(t)||LP(n;L2(g2)) < Cp • 

for all p > and k > Iq, where Cp is a constant that depends on ||Xo||^max{p,2)(f^.j^2(§2)), p, C, 
and a. 

Proof The result follows for p < 2 with Lemma 7.1 and with Holder's inequality. So we 
assume that p > 2 from here on. Let t £T and = to < ■ ■ ■ < tn = t he a partition of [0, t] 
for some n G N. We first observe that 



e=K+i m=-e 



LP(n;L2(S2)) 
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oo n—1 

+ 11 y ye-wEr=7+i/^»^,(_^-) 

II ^ ^ ^ ^-^^ LP(n;L2(S2)) 



Similarly to the proof of Lemma 7.1, the second term is equal to the norm of the approxima- 
tion error of an isotropic Gaussian random field with angular power spectrum {A^aJ^^ , £ G No) 
which satisfies by Lemma 5.3 that 

oo n—1 



l=K+l 3=0 



\a — 2 a—lJ 

Furthermore the first term satisfies similarly to the proof of Lemma 7.1 that 

oo i 

which converges faster than any polynomial and therefore can be bounded by So 
combining these two estimates, we finish the proof. □ 

Corollary 7.3. Let t G T and = to < ■ ■ ■ < tn = t be a discrete time partition for 
n G N, which yields a recursive representation of the solution X of Equation (2) . Furthermore 
assume that there exist £o £ a > 0, and a constant C such that the angular power spectrum 
{Ai,£ G No) satisfies < C ■ 1"°^ for all i > Iq. Then the error of the approximate solution 
X'^ is bounded uniformly in time, independently of the time discretization, and asymptotically 
in n by 

\\X{t)-X-{t)\\L^{S^)<^^-^ 

for all ^ < a/2. 

Proof. The proof is similar to the one for isotropic Gaussian random fields in Corollary 5.4 
bTit for completeness we include it here. Let j5 < a/2, then the Chebyshev inequality and 
Lemma 7.2 imply that 

F{\\X{t) - X«(t)|U.(s.) > n-^) < Kf^PE{\\X{t) - X-{t)\\l,^^,^) < . 
For all ^? > {a/2 — the series 

oo 

K=l 

converges and therefore the Borel-Cantelli lemma implies the claim. □ 



In Figure 6, we show the corresponding error plots to Figure 5 but instead of a Monte 

Carlo simulation of the approximate ^^(11: L^(§^)) error we plotted the error of one path of 
the stochastic heat equation. The convergence results coincide with the theoretical results in 
Corollary 7.3. 
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Tiber of series elements k 



number of series eiements k 



(a) Angular power spectrum with parameter a — 1. (b) Angular power spectrum with parameter a = 3. 




number of series elements k 



(c) Angular power spectrum with parameter a = 5. 



Figure 6. Error of approximation of a path of the stochastic heat equation 
different angular power spectrum of the Q- Wiener process. 



Appendix A. Interpolation spaces 

In this appendix we give a more detailed introduction to interpolation spaces than in 
Section 3 and show uniqueness of the spaces, i.e., that they are independent of the chosen 
interpolation couple. 

We consider the sequence of spaces 1, l),n £ No) that was introduced in Section 3 and 

start now with the definition of fractional order spaces by the real method of interpolation (see, 
e.g., [23, Chap. 1]). We observe that for any two integers k,n gNq with < k < n the pair 

is an interpolation couple with ^"(-1,1) C V^{-1,1) C L'^{-1,1). 
For integers k,m,n G Nq with 0<k<m<nso that < 6 := {m — k) /{n — A;) < 1, we may 
therefore define the intermediate space at "fine-index" q € [1,+cxd] 

s™f'"^(-i,i) = (y'=(-i,i),y"(-i,i)) 
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by the real method of interpolation as is introduced in [23, Chap. 1]. Then, these spaces are 



equipped with the usual norms || • || pm,(fc,n). given by 



"11^2- (-1,1) - \ 



poo 

J t-'^'i\K{t,u)\ij] , for l<q<+oo, 



supooi K(t,u), for ^^ = +00, 

where the if-functional is defined by 

^(*' = J^l^ {\Mv^{-i,i) + i|kl|y"(-i,i)) 

for t > 0. We observe that in particular the pair of spaces {V"'{—l,l),V"'~^^{—l,l)) is an 
interpolation couple for every n G No- Therefore, with n G No and for 1 < g < 00, we 
may extend the family {B^^^^^'^''\—l, l))o<fc<m<n,ge[i,+oo] of exact interpolation spaces also to 
noninteger numbers s = n + 9, 9 e (0,1), via 

Let us from here on simplify the notation and denote V'^{—1, 1) by and Bl^'^^'^\—1^ 1) 

by g*"^'"'' • Our next proposition states that for g = 2 and m G N, the Besov spaces B^2^'^'' 
are equal to for any choice k < m < n. 

Proposition A.l. Let m G N 6e given. For any k,n eNq with < k < m < n, it holds that 
r>22 — V ■ 

Proof. This result is classical (sec, e.g., [23], [18] or [20, Chap. 6.5] and the references there). 
We present the detailed argument here for completeness. 

By Proposition 3.2 we already know that the norm in is equivalent to weighted square 
summability of the coefficients of the Fourier-Legendre expansion. So it is sufficient to show 
the equivalence of the i?™2^'^'"^-norm and the convergence of the sum for all < A; < m < n. 

Therefore wc choose any fc, n G Nq with < k < m < n and u G 1, 1). Then u admits 
the Fourier-Lcgcndre expansion 

as has been seen above. Consider now u G V"^ U 3^2'^'"'^ C V'^. We split u into the sum 
V + 'w with V eV'^ and w eV" and the series expansions 

x2£ + l , ^2^ + 1 

V = 2_^{ue - we) — - — Pe and w = wi — - — Pe. 

e=o £=0 
Then Proposition 3.2 implies that 

K{t,ufc^ inf i\\v\\lk +t'^\\w\\ln) 
2^ + 1 



u=v-\-w 
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We observe further that the infimum over allu = v + w is equal to the infimum over all square 
summable sequences (lu^)^gNo ^ ^^(No), i-e., 

K{t,uf c^i inf V — - — Ge{ue,wf,t,k,n), 

where 

Ge{a, d:, t, k, n) := {a - df{l + f") + t^d\l + f"") 

is with respect to d G M a quadratic polynomial with positive leading coefficient for all ^ G Nq. 
For £ G No, its minimum is attained at 

where 

1 + 

This implies that 

"2^+1 



2^+1 2/-I , /,2A:\.2 9krM) 



t=0 
oo 



- - + f'gkn{e) 

and leads with the definition of the norm and the theorem of Fubini-Tonelli to 



roo 

e—r, 2 Jo 



- -+t^9kn{i) 



dt 



dt. 



where := {m — k)/{n — k) G (0, 1). To finish the proof it remains to show that 
by Proposition 3.2, i.e., we have to prove the equivalence 

i-oo ^1-29 

(1 + f'')9U^) / ^ 1 + ^ 1 + ^2((l-e)fc+en)_ 

Therefore let us split the integral first into 

POO ^1-26 rgkni^r^^'^ fl-2e roo ^1-29 

Jo i + t^gUt)'^^^ Jo 1 + t2fffc„(£) '^^ + 1 + f^gU^) 

and bound the two terms on the right hand side from below and from above by 



dt 
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and 

- — — / r'^~'^^dt< / ^ —dt 



This implies overall that 



1.1-261 



1 /"°° 



and moreover that 

/•oo ^1-26 

(l+^'")5fenW 1 + t2g^^(^) dtc^{l+ e'')gkn{^)' = (1 + + 

We observe that the function x^, p G (0, 1) is concave on M+ and satisfies (x + y)^ > 2*'~^(a;^ + 
yP), which implies finally that 

(1 + £2fc-)l-0(j ^ ^2n)e ^ (1 ^ ^2(l-0)fe-)(j ^ ^2en-) ^ ^ ^ ^2((l-0)fc+en) ^ j ^ ^2m_ 

This concludes the proof. □ 

Based on Proposition A.l, it is clear that one can use for every m G N in place of 5^2*^'^'"^ 
simply V^. Moreover, for fractional rj = n + 9 with n G No and some < < 1, we write 
also in place of ^2 2 • 
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